Sign change of the Griineisen parameter and magnetocaloric effect near quantum 

critical points 
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We consider the Griineisen parameter and the magnetocaloric effect near a pressure and magnetic 
field controlled quantum critical point, respectively. Generically, the Griineisen parameter (and 
the thermal expansion) displays a characteristic sign change close to the quantum-critical point 
signaling an accumulation of entropy. If the quantum critical point is the endpoint of a line of 
finite temperature phase transitions, T c oc [p c — p) , then we obtain for p < p c , (1) a characteristic 
increase r ~ T- 1/(vz) of the Griineisen parameter V for T > T c , (2) a sign change in the Ginzburg 
regime of the classical transition, (3) possibly a peak at T c , (4) a second increase V ~ — J— 1 /^ 2 ) 
below T c for systems above the upper critical dimension and (5) a saturation of T oc l/(p c — p)- 
We argue that due to the characteristic divergencies and sign changes the thermal expansion, the 
Griineisen parameter and magnetocaloric effect are excellent tools to detect and identify putative 
quantum critical points. 

PACS numbers: 71.27. +a 71.10.Hf 73.43.Nq 75.30.Sg 



I. INTRODUCTION 

The competition between two different ground states 
at a quantum phase transition leads to novel behavior 
in thermodynamics as well as in transport. A prominent 
example are the heavy fermion compounds whose non- 
Fermi liquid behavior is attributed to the presence of a 
quantum critical point (QCP) associated with a magnetic 
instability. 

We recently pointed oulii that the Griineisen param- 
eter, r, is an important tool to identify and classify a 
QCP since it necessarily diverges near a pressure-driven 
quantum phase transition with characteristic exponents 2 . 
In this article we focus on another aspect, i.e., the sign 
change of the Griineisen parameter. We will argue that 
generically the sign of Y (and therefore the thermal ex- 
pansion) changes as entropy is accumulated near a quan- 
tum critical point. The sign change in combination with 
the divergence leads to strong signatures of the Griineisen 
parameter near a QCP. 

A quantum phase transition (QPT) occurs at zero tem- 
perature upon tuning an external parameter like doping, 
pressure, electric field, magnetic field, etc. to a critical 
value. The underlying QPT manifests itself at finite tem- 
peratures in an unusual sensitivity of thermodynamics on 
these tuning parameters. In the following we will focus on 
QCPs which are tuned with pressure, p and/or magnetic 
field, H. Generalizations are straightforward. At T = 
the distance to the QCP is determined by the control 
parameter which depends on pressure and/or magnetic 
field, r = r(p,H). Near the QCP the control parame- 
ter can be linearized around the critical pressure or the 
critical magnetic field, 



r(p, H) = (p- Pc )/p = {H- H c )/H 



(1) 



where po and Hq are a constant pressure and magnetic 
field scale, respectively. Generically, the critical magnetic 



field will depend smoothly on pressure, H c — H c {p), and 
vice versa. Note that a linearization in magnetic field is 
only possible if the critical field is large, H/H c <C 1; in 
particular, this is not fulfilled in the case of a zero-field 
QCP. 

The critical contribution to the free energy density is 
a function of this control parameter and temperature, 
/ = f(r,T). The sensitivity on the tuning fields is ther- 
modynamically measured by the derivatives of the free 
energy density with respect to r. For a pressure-tuned 
QCP an example is provided by the thermal expansion, 
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where V m is the molar volume. With Eq. follows 
that the thermal expansion is directly proportional to 
the mixed derivative d 2 f/(drdT). When the QPT can 
be tuned by magnetic field the same derivative can be 
accessed by the thermodynamic quantity (dM / dT) h ■ If 
the QCP is sensitive to both, p and H, the critical part of 
the thermal expansion and of the temperature derivative 
of the magnetization are expected to be proportional to 
each other near the QCP, 



dHr 



(dM / dT) h dp 



(3) 



Their ratio gives the dependence of the critical magnetic 
field on pressure. Similar relationships hold for the mag- 
netostriction, compressibility and differential susceptibil- 
ity which all yield d 2 f/dr 2 . 

The main quantity of our interest is the Griineisen pa- 
rameter which is measured by the ratio of thermal ex- 
pansion and molar specific heat C p 
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FIG. 1: Different regimes in the phase diagram of a quantum 
phase transition. The dotted lines correspond to crossovers 
between the low-T and the quantum critical regime, T ~ \r\ vz . 
The control parameter might be sensitive to pressure and/or 
magnetic field. The solid line shows a generic isentrope along 
which the entropy is constant, dS — 0. 



Note that wc define the Griineisen parameter with the 
specific heat at constant pressure^ and not at con- 
stant volume (as often used in the literature). It is the 
specific heat at constant pressure which is measured in 
experiments on pressure-controlled QCPs. The corre- 
sponding quantity for magnetic field tuning is given by 
T H = - {dM/dT) H /C H , where C H = T^\ R . It will be- 
come clear below that Th describes the magnetocaloric 
effect. 

Let us shortly recapitulate the main results of Ref. 
The main observation is that T diverges at a quantum- 
critical point, while it is finite in all non-critical systems 
or close to generic classical critical points. This can most 
easily be seen in cases where scaling applies (i.e. for 
systems below the upper critical dimension, d + z < 4) 
and where the qualitative behavior of T can be extracted 
from the scaling form of the free energy 



f(r,T) = b-^ d+ ^f(rb x ' v ,Tb z ) . 



(5) 



where b is an arbitrary scaling parameter, d is the di- 
mensionality and v and z are the correlation length and 
dynamical critical exponent, respectively. As can be read 
off directly from Eq. Q , T scales like 1 jr or equivalently, 



dim[T cr ] = — dim[r] = —1/v. 
Accordingly, one obtains directly from JSJ) 

r cr oc — — ; — . 



(6) 



(7) 



in the quantum-critical regime, i.e. for T \r\ uz (see 
Fig. On the other hand, in the two low-temperature 
regimes on the right and left hand side of the QCP in 



Fig-d the Griineisen parameter diverges with the inverse 
of the control parameter r oc p — p c , 
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1 



V m (p-Pc) 



(8) 



Surprisingly, due to the third law of thermodynamics, 
i.e. by assuming a vanishing residual entropy at zero tem- 
perature, it is possible to determine even the prefactor G r 
of the divergence from a scaling analysis. It is given by 
a simple combination of critical exponents 



G T — 



y z-d 



yf 



(9) 



where the exponents i/q and y are associated with the 
low-temperature behavior of the specific heat, C p ~ T y o , 
on the right and left hand side of the QCP, respectively. 
As was shown in Ref. 01 these results might even hold (up 
to possible logarithmic corrections) in situations where 
the simple scaling Ansatz (JjjJ fails, i.e. for systems above 
the upper critical dimension. 

Equation © implies not only a divergence of T but 
also a sign change (assuming that G r has the same sign 
on both sides of the QCP)! Obviously the question arises 
where and how this drastic sign change takes place in 
the finite-temperature phase diagram. This will be one 
of the main topics discussed in this paper. 

The following section will discuss the sign changes us- 
ing qualitative arguments. Sec. IIIII investigates quan- 
tum critical points where there is no phase transition 
at finite temperature and briefly discusses experiments 
close to metamagnetic quantum phase transitions. In 
Sec. II VI we study how thermal expansion, Griineisen pa- 
rameter and magnetocaloric effects are influenced by a 
phase-transition at finite T in proximity to a QCP. An 
overview of our main results is given in Sec. IVl 



II. SIGN OF THE GRUNEISEN PARAMETER 

In order to obtain insight into the meaning of the sign 
of the Griineisen parameter it proves useful to consider a 
line of constant entropy within the pressure-temperature 
plane (p,T), 
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(10) 



Using the definition of the thermal expansion and the 
specific heat we obtain for T, 



1 dT 



V m T dp 



(11) 



The Griineisen parameter measures the variation of tem- 
perature upon pressure changes under constant entropy 
conditions. The Griineisen parameter thus corresponds 
to a pressure-caloric effect. As already alluded to, for 
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a QPT that can be driven by magnetic field the quan- 
tity analogous to the Griineisen parameter is the magne- 
tocaloric effect 



H 



(dM/dT) H _ 1 dT 
~ T dH 



C. 



H 



(12) 



where Ch is the specific heat at constant H . Experimen- 
tally, the quantities T and T h can be directly accessed by 
measuring the change in temperature at constant entropy 
upon pressure and magnetic field variations, respectively. 
In mathematical terms both yield the slope of the con- 
stant entropy curves, i.e. isentropes in the phase diagram. 

How do the isentropes look like near a quantum phase 
transition? We expect that we have an accumulation of 
entropy near the quantum critical point since directly at 
the QCP the system is frustrated of two different possi- 
ble ground states. From this expectation follows that the 
isentropes are tilted towards the QCP with a minimum 
in its vicinity, see Fig. ^ The minima of the isentropes 
indicate how the entropy accumulates around the QCP 
and sit at the positions where the system is maximally 
undecided which ground state to choose. According to 
Eq. (jl 1|> , the Griineisen parameter is proportional to the 
slope of the isentropes, i.e. it has a different sign on each 
side of the QCP. It thus follows from a generic entropy 
distribution around the QCP that the Griineisen param- 
eter will change sign in its vicinity. The sign change oc- 
curs at the location of the isentrope minima and therefore 
tracks the accumulation of entropy in the phase diagram. 
In this way, the Griineisen parameter maps out the en- 
tropy landscape near a quantum critical point. 

Since the specific heat is always positive the sign 
change of the Griineisen parameter coincides with that 
of the thermal expansion. As the thermal expansion is 
given by the negative derivative of the entropy with re- 
spect to pressure, a oc —dS/dp\T, a negative thermal 
expansion is present whenever the entropy increases as 
a function of pressure. This happens naturally in the 
vicinity of a pressure-controlled QPT. The sign change 
of the thermal expansion and hence of the Griineisen pa- 
rameter occurs at a pressure value where the entropy 
reaches a maximum. Similarly, for magnetic field tuning 
the magnetocaloric effect and the quantity (dM / 8T)h 
will change its sign at the accumulation point of entropy. 

As mentioned above, already the scaling result for the 
Griineisen parameter in the low-temperature regime, JSJl, 
suggested a sign change of T provided that the prefactor 
G r has the same sign in both low-T regimes. That this 
is the case is again ensured by entropic constraints. The 
sign of G r is determined by the relative size of the expo- 
nents y^ and dj z. These exponents determine the behav- 
ior of the specific heat or, alternatively, the entropy, S, 
in the low-T regime, S ~ T y o , and the quantum critical 
regime 1 , S ~ T d / Z , see Fig.^ Using again the argument, 
that the competition between two different ground states 
leads to an enhanced entropy close to the QCP, we ex- 
pect on physical grounds that the entropy as a function 
of T decreases in the low-T regime at least as fast as in 



the quantum critical regime. This implies that d/z < y^ 
and finally 



G r < 



(13) 



in both low-temperature regimes. There are examples of 
critical theories with exponents d/z — yo, so that G r = 0, 
e.g. an insulating Heisenberg antifcrromagnct with z = 1 
and yo — d on the ordered side of the phase diagram or 
an itinerant antiferromagnet with d = z = 2 and yo = I. 
In the latter example, logarithmic corrections ensure that 
the entropy at the critical point is higher 1 . 

In the low-temperature regime we can determine the 
critical isentropes explicitly by using the scaling result 
Eq. ©. We obtain for T <C \r\ uz 



T(r) 



(14) 



The isentrope behaves as a powerlaw near the QCP with 
an exponent given by G r , Eq. @. A minimum of the 
isentrope directly follows if G r is negative^. 

In the following we will give examples of two possible 
scenarios. If the entropy landscape at finite temperature 
is only determined by the zero-temperature transition we 
expect that the minima of the isentropes are located ap- 
proximately above the QCP, i.e. at r ~ 0. In this sce- 
nario, considered in Sec. IIIII the Griineisen parameter 
and the magnetocaloric effect change their sign near the 
critical pressure p c and critical field H c . If the QCP is 
however an endpoint of a line of classical finite tempera- 
ture transitions we expect that the entropy landscape is 
distorted with the minima in the vicinity of the critical 
temperature as sketched in Fig. We will show that 
the sign change of the Griineisen parameter then occurs 
within the Ginzburg region of the symmetric phase. This 
scenario is considered in Sec. IIVI 



III. SIGN CHANGE OF V NEAR A QCP 
WITHOUT ORDER AT FINITE T 

A. Ising chain in a transverse field 

A simple example of a QPT where the entropy land- 
scape is solely determined by the underlying QCP is pro- 
vided by the model of an Ising chain in a transverse field. 
The Ising chain shows a QPT as a function of transverse 
magnetic field H between a magnetic and a paramag- 
netic ground state. We are interested in the behavior 
of the magnetocaloric effect (|12f) near the critical field 
H c , which was also considered in Ref . @. The continuum 
theory describing the transition is given by (Majorana-) 
fermions with a relativistic spectrum^ 



(15) 



where r oc H — H c . The important exponents of the 
critical theory are z = d = v = 1. Furthermore, the 
spectrum is gapped away from the QCP which leads to a 
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FIG. 2: Sketch of the isentropes near the QCP of the Ising 
chain in magnetic field with r oc H — H c . The entropy ac- 
cumulates above the QCP leading to a sign change of the 
magnetocaloric effect 11211 at r = 0. 



specific heat that decays exponentially with temperature. 
The exponent i/q appearing in (j5J) can effectively be set 
to infinity on both sides of the QCP, so that the prefactor 
in both low-temperature regimes simplifies to 



-1 



(16) 



The thermodynamic quantities can be computed from 
the free energy density 



f(r,T) = -T 



dk 
2^ 



log 



2 cosh 



2T 



(17) 



The isentropes near the critical field H c are sketched in 
Fig. |21 Since the control parameter enters the free energy 
only quadratically the entropy landscape is symmetric 
with respect to the reflection r — » — r. This symme- 
try is rooted in the self-duality of the theory describing 
the QPT— . This has the consequence that in the case of 
the Ising chain the sign change of the magnetocaloric ef- 
fect occurs directly at r — 0. The self-duality symmetry 
thus causes the prefactor of the divergence Q to vanish 
in the quantum critical regime. The resulting magne- 
tocaloric effect is shown in Fig.[3]for a temperature sweep 
at constant magnetic field and vice versa. In the low- 
temperature regimes the magnetocaloric effect diverges 
with the universal prefactor i|16[) as is shown in the in- 
set of Fig. [3^. The developing divergence combined with 
a sign change leads to a strong signature of the magne- 
tocaloric effect in the field sweep shown in Fig. EJd. 



B. Experiments on metamagnetic quantum 
criticality 

An especially interesting example of a QCP without 
an associated finite-T phase transition are the so-called 
"metamagnetic quantum critical endpoints"— . The QCP 
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FIG. 3: Magnetocaloric effect 1121 of the Ising chain near 
the critical field H c as a function of the control parameter 
r — (H — He) /Ho for (a) a temperature sweep at constant 
magnetic field and (b) a magnetic field sweep at constant 
temperature. As shown in the inset, the saturation value in 
the low-T regime has the universal prefactor, —G r = vx = l. 



is here understood to be the endpoint, (H C ,T*), of a 
line of first order transitions where the temperature T* 
is tuned to zero. As the control parameter (i.e. the mag- 
netic field) couples here linearly to the order parameter 9 , 
the quantum-critical properties of such systems differ in 
some aspects from the examples considered in this paper 
- for a detailed discussion we refer to Ref. The no- 
tion of metamagnetic quantum criticality 9 was originally 
motivated by experimentsii*^ on Sr3Ru207- The meta- 
magnetic transition is sensitive to pressure and magnetic 
field variations suggesting that both p and H can be used 
as tuning parameters, r = r(p, H), and relations such as 
Eq. J2J) are expected to hold. In fact, the differential sus- 
ceptibility and magnetostriction nicely track each other 13 
near the critical field, H c s» 7.9T, confirming that pres- 
sure and magnetic field variations probe indeed the same 
thermodynamic information. The thermal expansion 10 
shows a change of sign near the critical field indicating 
the accumulation of entropy in the [H, T) plane above 
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H = H C . 

A metamagnetic anomaly of a different type but with 
qualitative similar signatures is also observed in the 
heavy- fermion compound CeRu2Si^ai£. The thermal 
expansion as a function of H again shows a pronounced 
sign change near the metamagnetic field H m ps 7.8T sug- 
gesting the vicinity to a QCP. According to Eq. (JSJ, the 
zero-temperature limit of the Gruneisen parameter is ex- 
pected to diverge as T CT (T = 0) oc 1 / (H — H c ) which also 
seems to be compatible with experimental observations, 
see Fig. I 6 of Ref. Q 



IV. SIGN CHANGE OF V NEAR A QCP WITH 
ORDER AT FINITE T 



It is a common situation that the symmetry-broken 
phase of the QPT extends to finite temperature, T. The 
QCP is then in fact the endpoint of a line of classical tran- 
sitions along which derivatives of the free energy show a 
singular behavior at a finite critical temperature T c . En- 
tropy is expected to accumulate near the phase boundary 
and the entropy landscape is correspondingly distorted as 
shown in Fig. 

As a consequence, the Gruneisen parameter is ex- 
pected to change its sign near the critical temperature 
T c . This is observed for example in the heavy fermion 
compound CeCug-z Au^r^ and in the spin-gap compound 
TlCuCls 22 . An especially nice example is the supercon- 
ducting heavy fermion system URu2Si2- The thermal 
expansion shows a pronounced jump at the supercon- 
ducting transition at T c = 1.18K accompanied with a 
change of sign 17 suggesting the vicinity of an associated 
QCP. The superconducting condensate of URu2Si2 forms 
in fact within a so-far unidentified "hidden order" phase 
which can be suppressed by an applied magnetic field of 
H c « 35. 9T 18 . Moreover, an additional reentrant phase 
is located in a magnetic field range of H = 36 — 39T. The 
magnetocaloric effect has been measured^ and the isen- 
tropes around the second order transition at H w 36T 
have a similar shape as in Fig. El with minima near the 
critical temperature, T c . Another example is the heavy- 
fermion alloy U(Pt,Pd)3 where a Gruneisen parameter 
inversion was also observed 19 . 

As an illustration of such a scenario we will consider 
the quantum phase transition of the dilute Bose gas. This 
example will capture many features which we believe are 
generic for a QCP with order at finite temperature. We 
assume that the chemical potential is sensitive to pres- 
sure so that a quantum phase transition can be induced 
by tuning the pressure to a critical value p c . More- 
over, for dimensions d > 2 a finite temperature transi- 
tion is present which distorts the entropy landscape in 
a manner as sketched in Fig. 0] The dilute Bose gas is 
relevant for the field-driven QPT in the spin-gap com- 
pounds which is interpreted as the Bose condensation of 




control parameter r 

FIG. 4: Sketch of a phase diagram with a symmetry-broken 
phase extended to finite temperature; here the QCP is the 
endpoint of a line of classical second order transitions. The 
minima of the isentropes, which identify the position of the 
sign change of the Gruneisen parameter, are located within 
the Ginzburg regime in the symmetric phase. Furthermore, in 
the Ginzburg regime the isentrope tend to nestle to the phase 
boundary as explained in the main text. For later reference 
the different regimes were labeled as: I/I' low-temperature 
regime, T <C \r\ uz , and II/IP quantum critical regime, T S> 
|r|" z , in the symmetric/ordered phase, and III the Ginzburg 
regime that houses the phase boundary, T c ~ (— r)*. 



The action of the dilute Bose gas has the form 7 



SI 



dr / d d r 



4>*{t, r) 



d 



V 2 - mo 



-||0(-r,r)l 4 



(r,r) 
(18) 



We will limit ourselves to a discussion of this model above 
the upper critical dimension of the QPT, d + z > 4 and 
therefore d > 2, where the correlation length and the dy- 
namical exponent of the QPT are given by their Gaussian 
values, 



I 

" = 2> 



(19) 



First let us consider the theory on the mean-field level, 
however, including the one-loop correction to the chemi- 
cal potential. The Landau potential takes the form 



(20) 



where the effective mass R is temperature dependent as 
it is renormalized by the critical fluctuations, 



R = 



-Mo 



2u 



d d k 1 , /k 2 
- — —r — coth — 
(2?r) d 2 V 2T 



magnons 



20,21,22,23.21.2") 



r - r cr (T) . 

(21) 
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In the last equation we introduced the control parameter 
r given in terms of the renormalized zero-temperature 
chemical potential, r = —fi = — /io + u J S^fa ■ If 
the QCP is sensitive to pressure changes, r can be lin- 
earized in the applied pressure, r = (p — p c )/po (assum- 
ing a finite magnetic field in the case of critical spin-gap 
compounds^), see Eq. JTJ. The phase transition occurs 
when the mass vanishes, R = 0. The phase boundary in 
the plane (r, T) is thus given by£ 



r C r(T) = 



((d/2) 

2 rf-l 7r rf/2 



(22) 



where the exponent reads 1/tp = (d + z — 2)/z = d/2. 
Note that the temperature dependence of the phase 
boundary, T c oc (p c — py, is due to the dangerously ir- 
relevant quartic coupling u. 

Above the upper critical dimension, d + z > 4, irrele- 
vant couplings will in general invalidate the simple scaling 
form of the critical free energy density given in Eq. JSJ. In 
the present example the dependence on the quartic cou- 
pling, u, must be incorporated into a generalized scaling 
form 



f(r,T,u) = b-( d+z '>f(rb 1/u ,Tb z ,ub d+z - 4 ) 



(23) 



The dependence on u may spoil the predictions for the 
Griineisen parameter drawn from Eq. This is for 

example the case within the wedge of the Ginzburg region 
indicated in Fig. 0] which contains the phase boundary. 
Bearing this in mind we now go beyond the mean field 
treatment by including the contribution of fluctuations 
in the free energy density. 



A. Gaussian approximation 

Depending on the regime in the phase-diagram plane, 
Fig. the physics of the dilute Bose gas is determined 
by different fixed points 2 ^. 

In the low-temperature regime I' within the symmetry- 
broken phase, T < — r, the physics of the Griineisen 
parameter will be controlled by the Goldstone modes; 
an expansion around the Gaussian theory in the quar- 
tic coupling u is plagued with IR divergencies in this 
regime indicating the crossover to the stable Goldstone 
fixed pointp2£i2L2&. This regime is conveniently treated 
by decomposing the fluctuations into massive amplitude- 
and gapless phase modes. The thermodynamics will be 
determined by the phase modes leading to a specific heat 
which decreases in temperature as C p ~ (T/y/^r) v ° with 
t/o = d. The prefactor, G r , of the Griineisen divergence 
© is here given by 



-G r \ 



u(d - zy ) 
Vo 



(24) 



On the other hand, the thermodynamics in the low- 
temperature regime I within the symmetric phase, T < r, 



is determined by the Gaussian fluctuations around the 
T = theory. The specific heat decreases exponentially 
with temperature, C p ~ r d / 2 (r/T) 1 / 2 e _r / T , reflecting 
the gap in the fluctuation spectrum. Correspondingly, 
in this regime the prefactor G r is given by 



-G r \ 



1 



(25) 



The most interesting regime for our discussion is the 
quantum critical regime, T » \r\. Here thermodynam- 
ics is almost always (in regions indicated as II and IF) 
dominated by ideal gas behavior with a specific heat, 
C p - T d ' 2 and thermal expansion a ~ 

T d/2-i_ 0nly 

sufficiently near the classical critical transition in regime 
III, thermodynamics will be controlled by the classical 
critical fixed point which for the dilute Bose gas belongs 
to the XY-universality class. This crossover occurs at the 
Ginzburg temperature when the Ginzburg parameter is 
of order one, 



g = uT\R\( d - 4 ^ 2 = <D(l) . 



(26) 



Note that the effective classical quartic coupling is given 
by u T . The description of the classical critical properties 
within this Ginzburg regime is beyond the simple approx- 
imation scheme employed here. It is within the Ginzburg 
regime where the thermal expansion and the Griineisen 
parameter change sign. However, as we will explain in 
detail below the sign change is not a property associated 
to classical criticality but is rather to be attributed to the 
underlying quantum phase transition. The sign change is 
a property of the quantum critical background on which 
the classical singularities develop. 

For its description we will evoke a Gaussian approx- 
imation which captures the correct thermodynamics in 
the quantum critical regime, T ^> \r\, except in the 
Ginzburg region where it will fail yielding singularities 
with wrong Gaussian exponents. In the quantum critical 
regime within the symmetric phase the Gaussian fluctu- 
ations determine the free energy density 



T 



d d k 



(27T) 



1 — exp 



R{r, T) + k 2 
T 



(27) 



For the following it will be crucial that here and in 
Eq. (|2*S|) the temperature dependent mass R(r, T) enters, 
see Eq. (|21|) . which measures the distance from the clas- 
sical transition. Note that in Eq. I|21(l we have set R = 
on the right-hand side, using only critical fluctuations 
when computing the mass renormalization. For our dis- 
cussion it is essential that this approximation is justified 
outside of the Ginzburg regime, as corrections are of or- 
der of the Ginzburg parameter, 5R/R = 0(Q). Outside 
the Ginzburg regime, we therefore obtain similar results 
as other approximation schemes like the self-consistent 
Hartree-Fock or the Popov approximation. Those ap- 
proximations have, however, the disadvantage that they 
wrongly predict a first-order transition22iSl within the 
Ginzburg regime, while our approximation describes a 
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second-order phase transition with Gaussian (and there- 
fore wrong) critical exponents. 

In the symmetry-broken phase the field fluctuates 
around the solution of the mean- field potential (3Jj). The 
condensate then attains the finite value, \cf>\ 2 = —R/u, 
and contributes to the free energy. The fluctuations 
around the finite condensate have the Bogoliubov spec- 
trum, 



/- = 



R(R-2r) 
2u 
d d k 
(2n) d 



(28) 



T 



log 



1 — exp 



Vk 2 (-2i? + k 2 ) 



Note that by virtue of (|2I \ the resulting entropy derived 
within this Gaussian approximation is indeed continuous 
at the critical temperature as is appropriate for a second 
order phase transition. 

The temperature dependence of the free energies, Q27JI 
and H28fl . is two-fold: there is an explicit T-dependence 
and an implicit dependence via the effective mass R — 
R(r, T) . In the following it will be useful to distinguish 
between a quantum critical and a classical critical contri- 
bution to thermodynamics. We will define the quantum 
critical contribution to be the one which derives from the 
explicit temperature dependence of the free energies. The 
classical contribution to thermodynamics will arise from 
the implicit T-dependence via the effective mass R(r, T) 
which is induced by the dangerously irrelevant quartic 
coupling u, see Eq. (122ft . This latter contribution will 
dominate thermodynamics near the classical transition 
since near T c the free energy is very sensitive to variation 
of R. 

To illustrate this point let us rewrite the free energy 
densities and lf2%|). 

f± = V±(r,R) +T {d+z ^ z T ± {RT- 1 ^ v ^), (29) 

where = and *_(r,ii) = (R(R— 2r))/{2u) is the 
contribution from the condensate. The part of the free 
energies that depends explicitly on temperature obeys 
scaling with z = \ jv = 2 and the scaling functions 



where K^ 1 = 2 d - 1 n d / 2 T{d/2) and 



J±(t,x) 



(30) 



u) + (t,x) = t + x, u-{t,x) = yjt(t - 2x) . (31) 

The Gaussian approximation yields an expression for the 
free energy density which conforms to the general scaling 
form (|2"5|l . The dependence on the quartic coupling, u, 
however appears only via the thermal renormalization of 
the mass R — R(r,T). Apart from this implicit temper- 
ature dependence induced by the dangerously irrelevant 
quartic coupling, u, the expression l|29|l resembles the 



quantum critical scaling form JSJ. Thermodynamic con- 
tributions that only derive from the explicit temperature 
dependence will therefore conform with the results ob- 
tained from Eq. (jSJ). In this sense it is appropriate to 
call these contributions quantum critical. The implicit 
temperature dependence via R = R(r, T) results in ad- 
ditional classical contributions to thermodynamics which 
are subleading except in the Ginzburg regime where the 
Gaussian approximation breaks down. The purpose of 
including the thermal renormalization of the mass is to 
ensure the correct threshold behavior at T c . The sign 
change of R(r, T) across the phase transition will reflect 
itself in a sign change of the Griineisen parameter. This 
sign change persists even outside the Ginzburg region as 
is explained below. 



B. Sign change of the Griineisen parameter 

We will show in the following that the sign change of 
the Griineisen parameter stems from the quantum critical 
contribution to the thermal expansion and is in fact a 
property of the functions T±. Consider the contribution 
to the thermal expansion at the critical temperature T c 
deriving from the explicit temperature dependence, 



8 2 

lim »„„„ /± 



*Q° P ~ V m po /'■•-•' ORDT 
1 



(32) 



V m p 



v(d + z) - 1 , , . 



Although the scaling function itself is continuous at the 
phase transition, J~+(0) = J--(0), its derivative is not. In 
our specific example we have 



(33) 



i.e. the quantum critical contribution to the thermal ex- 
pansion is discontinuous and changes sign at the phase 
transition. Note that the quantum critical contribution 
to the specific heat is smooth at T c . 

Making use of Eq. Q21f) . the associated anomaly in the 
thermal expansion, Achqcp = Qqcp 
lated to the derivative of the phase boundary, 



a QCP' can b e r e- 



Aoqcp 



1 1 3r cr (T) 



V m p u dT 



R=0 



V m plu [ dp 



dT r 



(34) 



The same anomaly follows also from the renormalized 
mean- field potential (|20|l . In this sense the anomaly 
Aojqcp and the resulting sudden sign change of the 
Griineisen parameter near T c can be interpreted as a fi- 
nite temperature manifestation of the mean-field char- 
acter of the underlying quantum phase transition. Al- 
though the sharp jump will be smeared out by the clas- 
sical critical fluctuations in spatial dimensions d < 4, 
the smearing is confined to the Ginzburg regime which 
is vanishingly small near the QCP. A pronounced jump 
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near T c in the thermal expansion and in the Griineisen 
parameter with an accompanying sign change will result. 

The smearing of the jump will also shift the exact po- 
sition of the sign change away from the critical temper- 
ature T c . In the following we will argue that this shift is 
towards the symmetric phase. 



C. Location of the sign change 

Where is the position of the sign change and hence the 
minima of the isentropes exactly located? Let us first 
give some general arguments. Let us consider the behav- 
ior of the entropy upon approaching the phase bound- 
ary from the ordered phase, see Fig. 0] The entropy 
attributed to the QCP increases when the phase bound- 
ary is approached from the ordered side by increasing 
the control parameter, dS/dr\r<T c > 0. When we en- 
ter the Ginzburg regime the change in entropy becomes 
dominated by the hnite temperature phase transition. 
The symmetric phase can be entered by either increasing 
control parameter r or temperature T. However, within 
the Ginzburg regime the tuning of r, i.e., pressure or 
temperature have the same effect since both parameters 
couple to the same relevant operator of the classical tran- 
sition. Since the entropy always increases as a function 
of temperature we also have dS/dr\x=T c > 0. The en- 
tropy as a function of r should therefore attain its max- 
imum, dS/dr = 0, above the critical temperature T c . 
As a consequence, it follows that the sign change of the 
thermal expansion and the Griineisen parameter should 
occur within the symmetric phase. 

We can obtain an explicit expression for the Griineisen 
parameter at the critical temperature, T c , when the spe- 
cific heat is sufficiently singular at the classical second 
order phase transition. The derivation follows standard 
arguments^. Near the finite temperature transition the 
singular part of the molar entropy can be written in the 
form 



Scl = S C h(T - T c (p),p) . 



(35) 



Near the critical temperature, the leading contribution 
to the thermal expansion will derive from the pressure 
dependence of the first argument. We thus obtain for 
the critical thermal expansion 



V m a r . 



dS, 



CL 



dp 



dS C L dT c 
dT dp 



(36) 

T c dp 



If the classical critical contribution is sufficiently singular 
such that the background contribution can be neglected 
the Griineisen parameter at T c is just given by the slope 
of the phase boundary, 



r(r = t c 



1 dT c 
V m T c dp 



V m (p c -p) 



(37) 



T=T c (p) 



In the last equation we made use of T c oc (p c — p) . The 
negative slope of the phase boundary i.e. the suppres- 
sion of T c with increasing pressure, r = (p — p c )/po, as 




-150, 



0.01 



0.02 0.03 0.04 
temperature T 



0.05 



0.06 



FIG. 5: Griineisen parameter of the dilute Bose gas as a 
function of temperature in dimension d = 3 as derived from 
the Gaussian approximation to the free energy, 12711 and 128> . 
The sign change occurs near the classical transition. The 
shaded area indicate the Ginzburg regime where the employed 
Gaussian approximation breaks down. The inset shows the 
temperature evolution of the specific heat coefficient and the 
thermal expansion. The value of the control parameter was 
chosen r = —0.001, temperature is shown in dimensionless 
units and the quartic coupling has been set to u = 1. 



depicted in Fig. 01 results in a negative Griineisen pa- 
rameter at the critical temperature, T(T = T c ) < 0. The 
Griineisen parameter thus has the same sign at the crit- 
ical temperature as in the ordered phase. We again find 
that the sign change must occur within the symmetric 
phase. 

In order to locate the sign change we consider the ther- 
mal expansion of the dilute Bose gas within the symmet- 
ric phase. The quantum critical contribution reads 



«QCP 



1 d 2 f + (R(r,T),T) 



V mPo dRdT 

2uV m p oT \ ) 

In addition to q;qcp, the implicit temperature depen- 
dence via the effective mass, R, yields a classical contri- 
bution 



1 d 2 f+(R(r,T),T)dR 



VmP0 



OR 2 



dT 



We obtain for their ratio 



> 



«CL 

aQCP 



2u 



d 2 f+ 
dR 2 



0(uTR^/ 2 ) . 



(39) 



(40) 



The classical part of the thermal expansion has a sign 
opposite to the part attributed to the QCP. The classi- 
cal part therefore reduces the contribution q;qcp upon 
approaching the phase boundary from the symmetric 
phase. This finally leads to the sought-after sign change. 
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The two contributions are of the same order when the 
Ginzburg criterion Il2(ill is fulfilled. The position of the 
sign change of the Griineisen parameter is hence located 
within the Ginzburg regime of the classical finite temper- 
ature transition. Sufficiently near the QCP the Ginzburg 
regime is vanishingly small so that the position of the sign 
change of T almost coincides with the critical tempera- 
ture T c (p). 

The temperature evolution of the specific heat, ther- 
mal expansion and the Griineisen parameter is shown in 
Fig. The control parameter has been chosen nega- 
tive such that the critical temperature is crossed during 
the temperature sweep. The temperature region shown 
in Fig. [3 corresponds to the quantum critical regime, 
|r| < T 1 /^), where the expressions (|27|) and l|28|) for the 
free energy are applicable. In this regime the Griineisen 
parameter is determined by temperature and obeys the 
scaling form J7J). Since we have vz = 1 it behaves as 
r oc X/T. The curve follows this behavior except within 
the small Ginzburg regime where the Griineisen param- 
eter strikingly changes its sign and forms a sharp peak 
at the critical temperature. Since within our Gaussian 
approximation the classical critical specific heat is diver- 
gent with the Gaussian exponent, acauss = 2 — d/2 > 0, 
the peak value is given by Eq. (|37fl . where for the di- 
lute Bose gas we have ^ = 2/d. In particular, at the 
phase transition the Griineisen parameter is now deter- 
mined by the distance to the QCP, p — p c . Since the 
phase boundary is located well inside the quantum criti- 
cal regime, T ^> \p — p c \, where T is usually determined 
by temperature, T oc 1/T, the Griineisen parameter is 
strongly enhanced at T c . The crossover to this enhanced 
value occurs within the narrow Ginzburg regime leading 
to a peak structure at T c . The peak of T at the criti- 
cal temperature manifests itself in an additional tilt of 
the isentropes within the Ginzburg regime as sketched in 
Fig. 0J Indeed, remembering that the Griineisen param- 
eter just measures the slope of the isentropes, Eq. ifTTJl. 
it follows from expression l|37l) that the isentrope locally 
follows the phase boundary at T c . 

The well-pronounced peak of the Griineisen parameter 
at T c is only expected for a sufficiently singular classical 
critical specific heat such that relation (|37[1 holds. The 
Gaussian approximation employed here overestimates the 
classical critical specific heat exponent of the dilute Bose 
gas. In fact, the specific heat exponent of the d = 3 
XY-universality class is negative, i.e. the specific heat 
is not divergent at T c . Nevertheless, we still expect it 
to dominate over the quantum critical background such 
that a narrow peak in T should evolve at T c . 

As can be seen in the inset of Fig. [31 near the QCP the 
anomaly in the specific heat, AC P , near T c is small in 
comparison with the one in the thermal expansion, Aa. 
This is expected since the large anomaly in the thermal 
expansion is attributed to the quantum critical point, 
Aa = Ackqcp, see Eq. (|34|) . whereas the anomaly of the 
specific heat originates only from the classical contribu- 
tion to thermodynamics due to the temperature depen- 



aV m 



P=Pc 



d + z > 4 




FIG. 6: Sketch of the thermal expansion coefficient, a/T, in 
different regimes of the phase diagram in Fig. H for a QPT 
above its upper critical dimension, d + z > 4. The ther- 
mal expansion changes sign near the critical temperature T c . 
The exponent a is the specific heat exponent of the classical 
transition and is determined by the spectrum of low-lying 
excitations in regime I'. 



dence induced in the effective mass, R, by the danger- 
ously irrelevant quartic coupling, u. 

This is also in agreement with the Ehrenfest relation^ 
which compares the size of the anomalies, i.e., the jumps 
in thermal expansion and specific heat at the classical 
transition that derive from the mean-field potential i|2U|l , 



Aa 1 dlog T c (p) 



4, 



A Co 



V„, 



dp 



V m {p-p c ) 



< 0. 



T=T c (p) 



(41) 



In the second equation we made again use of Eq. (|22ll . As 
the quantum critical point is approached the relative size 
of the anomalies is expected to diverge as oc l/(p — p c ) 
resulting in a dominant anomaly in the thermal expan- 
sion. Note that although it has the same functional form 
as Eq. (|37|l the Ehrenfest relation contains different in- 
formation. In particular, the Ehrenfest relation does not 
discuss the absolute size of the Griineisen parameter T. 

At low temperatures the Griineisen parameter will 
eventually saturate after crossing over into the low- 
temperature regime F (not shown) and converges to a 
value now given by Eq. (JSJ) . The behavior of the thermal 
expansion and the Griineisen parameter in the various 
regimes of the phase diagram are summarized in Fig. [B] 
andQ 



V. SUMMARY 

The Griineisen parameter, T, and the magnetocaloric 
effect change sign near generic quantum critical points. 
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FIG. 7: Sketch of the Griineisen parameter, F, (or equiva- 
lently of the magnetocaloric effect in the case of a field-driven 
transition) in different regimes of the phase diagram in Fig. ^ 
for a pressure-tuned QPT, r oc p — p c , above its upper critical 
dimension, d + z > 4, where v = 1/2 and = z/(d + z — 2); 
G r is defined in Eq. [5] Near the classical critical transition V 
changes its sign in a characteristic jump. The peak at the crit- 
ical temperature, T c oc (p c —p) , is present if the specific heat 
is sufficiently singular at T c , see Eq. 1371 . The intermediate 
regime IF vanishes in the limit d + z — > 4+. 

We showed that the position of the sign change indicates 
the accumulation of entropy in the phase diagram. Two 
scenarios have been distinguished: a QCP with and with- 
out a symmetry broken phase at finite temperature. For 
the latter, treated in Section ITTT1 the sign change is ex- 
pected to be located near the critical value of the control 
parameter, e.g., near the critical pressure p c . We dis- 
cussed two examples where such a scenario is realized: 
the Ising chain in a transverse field and metamagnetic 
quantum critical materials. 

An overview of our main results for the second sce- 
nario, where a symmetry-broken phase at finite T is 
present, is given in Figs. HH and [7| While explicit cal- 
culations have been performed for a QPT of the dilute 
Bose gas, we expect that all of our qualitative results are 
equally valid for other quantum phase transitions where 
the control parameter couples quadratically to the order 
parameter. Note that the exponents for T in Fig. were 
obtained by considering the ratio of the critical parts of 
thermal expansion and specific heat; for comparison with 
experiments a possible non-critical background contribu- 
tion might have to be subtracted, cf. discussion in Ref. 1. 



Most interesting is a situation where one investigates the 
behavior on the ordered side of the phase diagram close 
to the QCP. Upon lowering temperature one crosses ac- 
cording to Fig. |1| four different regimes. 

In regime II one observes the usual power-laws in 
thermal expansion and Griineisen parameter (or equiva- 
lently magnetocaloric effect) which are associated to the 
quantum-critical part of the phase diagram. The main re- 
sult of our paper is the characteristic jump of thermal ex- 
pansion and Griineisen parameter in the Ginzburg regime 
located slightly above the classical phase transition. The 
jump in r proportional to ^ > where VP charac- 

terizes the form of the phase boundary, T c oc (p c — p) , 
gets more and more pronounced upon approaching the 
QCP. If the classical specific heat is diverging, there is 
also a sharp maximum in T at T c with the universal value 
r(T c ) = iff/ {V m (p—p c )). In situations where the classical 
specific heat is not so singular, a peak can still occur but 
a lower maximal value is expected. In the ordered phase 
outside of the Ginzburg regime, i.e. in regime IF, the 
Griineisen parameter increases again as in II but with 
an opposite sign as the system is located on the other 
side of the phase transition. Finally, at lowest tempera- 
tures a saturation of T sets in to an universal value given 
by Eq. (JSJ. For a magnetic-field tuned quantum phase 
transition an analogous behavior is expected for the mag- 
netocaloric effect <|12[) . in cases where the QCP is below 
its upper critical dimension, one has iff — vz and region 
IF is absent. 

It is an interesting open question how the above results 
are modified for quantum phase transitions which involve 
different types of fluctuations possibly characterized by 
different time-scales as for example in the case of itinerant 
ferromagnetism^iiS. 

To summarize, the divergence of the Griineisen param- 
eter and magnetocaloric effecfei near a QCP in combina- 
tion with their sign change result in very strong signa- 
tures. They are thus important thermodynamic probes 
to detect and classify QCPs. 
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